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ABSTRACT 


Fracture 1S one of the dominant failure modes in 
structures subjected to external loads. Stress and deformation 
fields around the crack tip are important to understand the 
crack propogation and arrest. For a plate with a through-the- 
thickness crack and subjected to a bending load, there is 
crack closure on the compression side of the crack face. The 
present study investigates effects of crack closure on the 
stress and deformation fields on the tension side of the crack 
face. A three-dimensional finite element analysis is performed 
for laminated composite plates using both the line and surface 
crack closure models. For a composite whose longitudinal 
elastic modulus is much greater than the transverse modulus, 
line and surface closure models result in higher stresses near 
the crack tip in comparison to the no-closure solution. 
Hence, mno-closure solutions are nonconservative for the 
composite. Transverse shear is the major cause for the 


nonconservative solution. 
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I. INTRODUCTION 


Plates and shells are very commonly used as structural 
members in construction of ships and other structures. Due 
to their wide usage as structural members, much research has 
been performed to understand their failure mechanisms and 
modes under different loading conditions. Fracture is one of 
the most frequent failure modes in plates and shells. 

The first analysis on plates containing centrally 
located through-the-thickness cracks and subjected to 
bending loads was performed by Williams [Ref. 1] in 1951. In 
his analysis, Williams used the forth-order classical plate 
bending theory with Kirchhoff boundary conditions. He showed 
that there are some discrepencies between the stress 
distributions for extensional and bending problems. In 
William's bending study, surface shear stresses on crack 
faces didn't vanish. This physical violation is attributed 
to the fact that the Kirchoff bending theory can not satisfy 
the physical boundary conditions on the crack surface 
completely. 

This problem was solved by Knowles and Wang [Ref. 2]. 
They used Reissner/Mindlin's sixth order plate theory which 
allows three boundary conditions at edges. They solved the 


problem by employing a Singular integral equation. However, 


their solutions were limited to plates with very small 
thicknesses. Analyses including effects of the plate 
thickness on the bending stress distribution were performed 
independently by Hartranft and Sih [Ref. 3] and Wang [Ref. 
4}. 

All the theories mentioned above assume that the crack 
surface is stress-free and straight. In reality, for a plate 
under a bending load, there exists crack closure on the 
compression side and crack opening on the tension side of 
the crack surface. A photoelastic experimental study 
performed by Smith and Smith [Ref. 5] showed partial crack 
closure. A significant change of local stresses was observed 
as a result of the crack closure. 

A finite element analysis was first undertaken by Jones 
[Ref. 6] using Kirchhoff's bending theory to study effects 
of crack closure in plate bending. He compared the closure 
and no-closure solutions, and reported increases in crack 
opening displacements and intensity factors due to existing 
crack closure. ere Heming's study [Ref. ale the 
Mindlin/Reissner sixth order plate theory was employed, 
which is more realistic and allows satisfaction of three 
boundary conditions at the crack face. These results were 
contradictory to the results obtained by Jones ([Ref. 6]. 


Heming stated that there were significant decreases in the 


bending stresses and displacements on the tension side due 
emack closure. 

In those studies, the line closure model, which assumed 
the outermost points from the neutral axis on the 
compressive side of the crack surface to be in contact, was 
used. However, Smith and Smith's experiment [Ref. 5] showed 
that there was closure over an area of the crack surface, 
the size of which varied depending on the plate and crack 
dimensions. A three dimensional finite element analysis was 
performed by Alwar and Nambissan [Ref. 8] to investigate the 
closure region and effects on the stress distribution of 
this closure. A very crude finite element mesh was used in 
their analyses. 

Finally Kwon [Ref. 9] studied the crack closure effect 
for unidirectional composite plates using the line closure 
model and the Mindlin/Reissner plate theory. He found that, 
the crack closure on the compression side of a bent 
unidirectional composite plate increased the crack opening 
displacement and bending stress on the tensile side. This 
finding was in contradiction to that for an isotropic plate. 
Hence, it was concluded that results obtained without 
considering crack closure were nonconservative EGE 
unidirectional composite plates. The nonconservative 
solutions for unidirectional composite plates were caused by 


the effect of transverse shear. 


In the present study, effects of surface closure on 
stresses and displacements around the crack tip in laminated 
composite plates are investigated using three-dimensional 


finite element analyses. 


II. FINITE ELEMENT FORMULATION 


a; OVERVIEW 

A finite element program is developed to model a plate 
with centrally located through-the-thickness crack and 
subjected to a constant bending load. As shown in Figure 1, 
the cartesian coordinate system is employed in the model. 
Displacements in the x, y and z directions are represented 


by u, v, and w' respectively. 





Figure 1 Geometric representation of the plate model 


The plate is assumed to have simply supported boundary 
conditions along its edges. Due to symmetries along the x 
and z axes, only a quarter of the plate is analyzed to 
reduce the number of nodes needed in modeling the plate. 
Therefore, two symmetric boundary conditions and two simply 
Supported boundary conditions are required for the quarter 
plate. Three-dimensional parallel-piped elements are used in 
constructing the stiffness matrix. Each node has three 
degrees of freedom and, therefore, there are 24 dofs per 


element as shown in Figure 2. 





Figure 2 Nodal coordinates of a typical finite element 


Displacements of the plate are functions of x, y and z 


coordinates. That is, 


u=f\ (x,y, 2) 
v= f2(x,y, 2) (1) 
= = 3(x,y, 2) 


To interpolate the displacements within each finite element, 
Shape functions are utilized. Displacement can be expressed 


as, 


where subscript i represents the node number, and the 
numeric subscripts 1, 2, and 3 denote the displacements in 
aiemwee y, and 2 directions, ) respectively, and N. 15 “a 
mumeETOn of X, y and Zz. 

From the elasticity theory, strains are expressed as 


Given below, 


Cu aay _ ow 
ae a CV) > ey. Cz = 
i 
3) 
Cu. cv ee ce Pp ids ee i 
= a ax Cz aT Se Cyz = at & 


When written in the matrix form, 


(cligs = (Bieeoat one (4) 


where {e} is the strain vector of size 6 by 1, {u}©* is a 
nodal displacement vector of size 24 by 1, and [B] 1s a 
matrix containing derivatives of the shape functions. 

The three-dimensional generalized Hooke's law is 


defined as, 


{o} =[D]te} (5) 


where [D] is a 6 by 6 material property matrix, and {o} and 


{e} are stress and strain vectors. The stress and strain 


vectors are written as, 


{Oo} =[Oxx, Oyy, Ozz, Txy, Txz, esez| 


(6) 
{e} = [Exx, Cyys GZZz, Yxys Vez ‘eae 
The elastic strain energy can be written as, 
ae 
E=1| {e}"[D{ejav= + | (Us "(BI (DIB) Uyav (7) 
|” V 
In equation (7), only the terms in paranthesis have 


dependency on the coordinate system. The volume integral of 


these terms yields the stiffness matrix, 


[K]’ =| (2) [D][Blav (8) 
} 


The equation of static equilibrium for an element can 


be written as, 
ey = (K}]~{ 0} * (9) 


where {F}° is the vector of forces acting on the elemental 
nodes. 
B. CONSTITUTIVE EQUATION 

The material property matrix [D] is a symmetric matrix 
which includes elasticity moduli and Poisson's ratios of the 
material. In the present work, isotropic, orthotropic and 
laminated composite plates are examined. 

ane ISOTROPIC MATERIAL 

There are three independent material constants in 

the [D] matrix for an isotropic material. Elasticity moduli 
and Poisson's ratios are the same in every direction. For an 


1sotropic material matrix [D] is given as, 


E*(1-v) Fey Eeyp 
(l+ev)*(l-2*v) (ltu)*#(1-2*u) (1+u0)*(1-2#v) 0 0 0 
Fev E*(1l-v) Fev 0 
(1l+v)*(1-2*vu) (l+u)*#(1-2*vu) (14+v)#*(1-2*v) 
Etvu Etu E*v 0 
(1+0)*(1-2*v) (1+u)*(1-2*v)  (1+0)*(1-2*0) (10) 
0 0 0 GO 0 
0 0 0 0G 0 
0 0 0 OmO G 


Dis ORTHOTROPIC MATERIAL 
The number of independent material constants 
in the [D] matrix is six for an orthoropic material. The 
material properties have three principal axes. The material 


properties of the plate are considered as, E,/E,=5, 20, 40, 


and 100, G,./E,=0.6, G,,/E,=0.5, v,,=0.25, and 0 >j,=1/9. Imeem 
of the [D] matrix is called the compliance matrix and it is 


given in [Ref. 10] as shown below, 


W/E -02/E -031/E3 Q Q Q 
-V12/E | l/E> -32/E3 0 0 0 
—013/E —023/E 2 VWE3 0 0 Q (11) 
0 0 0 1/G23 Ss «0 0 
0 0 0 QO 1/G3, 0 
0 0 0 0 0 WGi~ 


where E denotes inplane elasticity modulus, G denotes shear 
modulus, and numeric subscripts 1, 2, and 3 4represent 


material principal axes. 


ILC) 


So LAMINATED COMPOSITE MATERIAL 

A laminated composite plate with [0/90]. is 
investigated. There are two [D] matrices in the model to 
represent the oriantations of 0 and 90 degrees. The matrix 
[D] for the 0 degree layer is the same as that for the 90 
degree layer except that material properties in the x and z 
directions (see Figure 1) are interchanged for the 90 degree 
layer. As a result, the compliance matrix for the 90 degree 


layer is given below, 


V/E3 —03/E 2 -013/E| Q Q Q 
—032/F3 l/E2 —V2/E| 0 0 0 
—03)/E3 —v>;/E 3 VE 0 Q 0 (12) 
0 0 0 1/G2, 0 0 
0 0 0 O 1/G;3 0 
0 0 0 0 OO G3, 


Cc STRAIN-DISPLACEMENT RELATION 
Linear shape functions are utilized in the present 


finite element formulation. That is, 
a-x b-y cz 

Ni = (2) (2) (=) 
x b-}' cz 

No =()*() «(= 


N= (3) *(F) 


ial 


Na =(3)*(%) «4 (13) 
Ns = t)* (7) «@ 
No =) * (2) * 


Na =H) *(Z) *® 


Ng = (3) * i) “(G) 
where a, b, and c are dimensions of a finite element as seen 
Tier VoUGeg 22 

Matrix [B], which relates strains to nodal 
displacements of an element is formed by taking partial 
derivatives of the shape functions. Hence, six strain 


components are expressed as shown below, 


oe) 


O 
Cx = » ali * Gin] 
8 
O 
Cyye= » ali * wp) 


g 
ex= DY LN; *Uj3] (14) 


IZ. 


8 8 
= 1S SIN, * Un] + ms £{N, * U2] 
a | =| 


8 8 
£IN,*UnjJ+ D FIN, * Us] 


me | =| 


Cx 


e= & sIN.*Unl+ £ ZIN.*Ua] 
j= 1 fed 


D. STIFFNESS MATRIX 
The element stiffness matrix [K]° is found by volume 


integration of the term ([B]’[D] [B]) as in the equation (7), 


where x varies from 0 to a, ( x € (o a})) 
Wawartes /Eromr0 to b, 4-y € [orb ) 
Zevarves fFroms0- C0 ¢€, ( zre .[o<c)) 


The global stiffness matrix 1s constructed by adding 
element stiffness matrices. In the present finite element 
mesh, the number of nodes in the x, y, and z directions are 
21, 9, and 21 respectively. It gives the halfbandwidth of 
636 and total degrees of fredom of 11907. A static pressure 
load is applied at the nodes ao the xz plane lying in the 


middle of the plate. 
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III. BOUNDARY CONDITIONS AND MODELING OF CRACK CLOSURE 


Selecting proper constraints imposed along the edges of 
a plate 1S very important. As mentioned in [Ref. 7], 
application Oe Ko jeloyoneie boundary condi tiem Oy 
Reisner/Mindlin boundary conditions results in contradictory 
solutions. Iie 1s more realistic ze) employ the 
Reissner/Mindlin theory, which describes three boundary 
conditions along the edges of a plate. Therefore, in the 
present models the Reissner/Mindlin BCs are utilized. 
Because of the symmetry involved, a quarter of the plate is 
modeled. The model has two symmetric and two simply 
supported boundary conditions. Symmetric BCs are satisfied 
by constraining displacements in the normal direction to the 
boundary surface. Simply supported boundary conditions, on 
the other hand are imposed by constraints in vertical 
displacements and tangential displacements at the boundary 
surface. 

Crack closure is a constraint on the crack face. aillew 
are three crack surface models which are used in analyses of 
plates with through-the-thickness cracks. The first model 
assumes that the plate has no constraint on the crack 


surface, as seen in Figure 3a. This 1s called the 


14 





Figure 3a No-Closure Model 





Figure 3b Line Closure Model 





Figure 3c Surface Closure Model 


Alig 


no-closure model. The crack surface is considered as 
stress-free and straigth 

The assumption of no-closure allows crack opening to 
occur on the tension side of the crack face. On the other 
side of the crack face, deformations of the same magnitudes 
occur in the opposite direction. Hence, overlapping or 
penetration takes place by employing the no-closure model. 
However, penetration is not possible in the actual case. In 
reality, there would be contact over a region on the crack 
surface. And this contact is called crack closure. 

We can model the closure region in two different ways. 
In the first model the outermost nodes on the compression 
Side with respect to the neutral axis are constrained to 
avoid penetration. This is called the line closure model and 
1s shown in Figure 3b. Alternatively, a region which 
contains some of the nodes between the outermost nodes and 
neutral axis on the compression side of the crack face can 
be constrained, as shown in Figure 3c. This second model is 
called the surface closure model. 

Plate bending theories permit only the line closure 
model, because the theories assume that plane sections 
remain plain before and after bending. The experimental 
Study by Smith and Smith [Ref. 5] shows that crack closure 
takes place over a portion of the crack face. Therefore, the 


surface closure model is the closest describtion of the 


1s 


crack closure. An iteration process 1S necessary to model 
Surface closure using a contact algorithm. Proper number of 
nodes are constrained to ensure that there 1S no penetration 
on the crack surface. 

When usSing the finite element method, it 1s possible to 
model the surface crack closure using 2-D plate bending 
elements or 3-D solid elements. If 2-D plate bending 
elements are used in modeling, a stack of plate elements are 
used through the plate thickness. As a result, the crack 
face can have a piecewise linear deformation and the closure 
portion becomes a rectangular shape. With 3-D solid 
elements, displacements at the crack face can have more 
general deformation. Hence, a general shape of surface crack 
closure can be modeled. Because of this advantage, we use 


3-D elements in the present study. 


aly 


IV. RESULTS AND DISCUSSIONS 


In the present study, cracked plates made of three 
kinds of materials are examined. These are isotropic, 
orthotropic and three layered composite plates with an 
orientation of [0/90/00]. Four different elasticity moduli 
ratios, E1/E2=5, 20, 40, and 100, are considered for 
composite and orthotropic models. Consequently, we analyzed 
nine models with different material properties. Three types 
of crack closure modelings, which are discussed in the 
previous chapter, are employed to model the crack surface. 
These are surface closure, line closure, and no-closure 
models. All three models are compared with each other for 
each cracked plate. 

All the stress and displacement data are computed on 
the plane containing the crack. A three dimensional 
representation of the crack surface for the no-closure model 
is shown in Figure 4. In this figure flat section is the 
uncracked portion of the crack plane. Since for the 
no-closure model there is no constraint on the crack 
surface, displacements with the same magnitudes and in the 
opposite directions are observed on the tension = and 


compression sides. In the surface closure model, due to the 
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crack closure, only displacements in the crack opening 
direction are allowed, as shown in Figure 5S. 

Effects of closure on the crack opening displacements 
can be explained by considering a compressive force due to 
closure as shown in Figure 6a, that acts on the compression 
Side of the crack surface. A statically equivalent force and 
moment, aS in Figure 6b, can be thought acting in the middle 
of the plane. While the force increases the crack opening 
displacements on the entire surface of the crack, the moment 
causes decreases on the tension side. For smaller elasticity 
ratios the effect of the moment iS greater. Hence, we 
observe decreases in crack opening displacements. For higher 
elasticity ratios, especially in the case of a composite 
plate, the effect of the force is higher. 

When the elasticity ratio is higher, shear deformation 
due to transverse shear becomes more important. As seen in 


Figure 6c the rotation caused by the moment is decreased due 


to the transverse shear effect. Here 8 is the rotation angle 


for the moment without considering transverse shear effect, 


and y 1s the rotation angle in the opposite direction due to 
the transverse shear. 

In the displacement’ plots, the horizontal axis 
represents the normalized distance from the crack tip and 
the vertical axis represents crack opening displacements on 


the tension side of the plate. These crack opening 


TENSION SIDE 


- Sa 


COMPRESSION SIDE 


Fisure sa Compressive Force due to Closure 


TENSION SIDE 


COMPRESSION SIDE 


Tigure 6b Statically Equivalent Force and Moment 


neutral axis 


displacements (CODs) are normalized by the maximum 
displacement of the corresponding no closure model. Since 
there 1S contact on the compression side for closure models, 
only displacements on the tension side are plotted. 

It is observed that the crack opening displacements for 
line and surface closure models are very close to each other 
for all the elasticity ratios. When low elasticity ratios 
are considered, crack closure decreases the crack opening 
displacements on the tension side all the way. This 1s due 
to the fact that the effect of the moment is higher for low 
elasticity ratios. Since there are more constraints for the 
surface closure model its displacements are smaller. Figure 
7 shows the crack opening displacements on the tension side 
of a laminated, cross-ply composite plate with the 
eWasticity ratio EL/E2=5. 

From Figure 8 it is seen that, when the elasticity 
ratio, EL/E2, is 20, crack closure reduces the displacements 
near the crack tip. On the other hand, i1t increases the 
displacements away from the crack tip. When the plots for 
E1/E2=5 and E1/E2=20 are compared, a slight increase near 
the crack tip is also observed. 

Further inreases in the crack opening displacements are 
obtained for elasticity ratio E1/E2=40, as seen in Figure 9. 
Near the crack tip, the displacements for closure cases are 


Still lower than that for the no cloSure caSe. 
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Figure 10 Crack Opening Displacements (CODs) on tne tension 
Side of a laminated, ross-ply composite sate with 
SuetoseaGliky ratio £,/E.=100. 


When E1/E2=100, closure increases the crack opening 
displacements all the way, as seen in Figure 10. This is 
because the transverse shear effect is high enough to 
overcome the effect of the moment also near the crack tip. 

Since the closure cases reduce the stresses on the 
compression side more, stress distributions on the tension 
Side are different from those on the compression side. 
Therefore, both tension and compression sides are examined. 
In the stress plots, the horizontal axis represents the 
normalized distance from the crack tip upto the corner of 
the plate along the uncracked portion of the crack plane. 
The vertical axis represents stresses that are normalized by 
the applied constant pressure . 

From Figures 11 and 12 it can be seen that for E1/E2=5, 
there exist large discrepencies in the stress levels between 
closure and no-closure cases. These differences occur both 
on the tension and the compression sides of the plate. On 
the other hand, the differences in stresses for surface 
closure and line closure cases are very small as seen in the 
crack opening displacements. This is especially true on the 
tension side. 

When higher elasticity ratios are considered, increases 
in the stresses near the crack tip are observed for closure 


cases compared to no-closure case, as shown in Figures 13 to 


18. When the E1/E2=100, higher stresses are obtained on the 


tension side for closure cases. 
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STRESSES ON THE TENSION SIDE 
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Figure ll Stress distribution on the tension Ssidemieuee 
laminated, cross ply composite plate with elasticity ratio 
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Stress distribution on the compression side of a 


cross-ply composite plate with elast 
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Figure 18 Stress distribution on the compression side of a 
laminated, cross-ply composite plate with elasticity ratio 
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V. CONCLUSIONS AND RECOMMENDATIONS 


When a laminated composite plate with a through the 
thickness crack has a large ratio between the longitudinal 
and transverse elastic moduli, crack closure increases the 
crack opening displacement and tensile bending stress around 
the crack tip compared to those obtained without considering 
crack closure. There is also a slight difference between 
surface and line closure solutions since the surface closure 
model has more constraints. Both crack closure solutions are 
nonconservative. AS a result, crack closure must _ be 
considered for most composite plates under consideration. 

An experimental study shows that the surface closure 
model is the closest description of crack closure. Since 
Similar displacements and stresses are, however, observed 
around the crack tip, line closure can be used to model 
crack closure for computational efficiency. 

In future studies experiments may be performed to 
validate the computer results. Effects of crack closure on 
composite plates containing delamination cracks need to be 


analyzed. 
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